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Abstract 

Hedgehog solitons are found in a chiral quark model with non-local interactions. 
The solitons are stable without the chiral-circle constraint for the meson fields, as 
was assumed in previous Nambu-Jona-Lasinio model with local interactions. 
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In recent years a lot of efforts have been undertaken to describe baryons as solitons of 
effective chiral models. Of particular interest are models which include the Dirac sea, 
since they allow for a unified description of mesons and baryons. So far solitons have only 
been obtained in the Nambu-Jona-Lasinio model with the proper-time or Pauli-Villars 
regularization of the Dirac sea [1-3] One problem encountered with the proper-time 
regularization, by far the most commonly used, is that solitons turn out to be unstable 
unless the sigma and pion fields are constrained to the chiral circle [5,6]. Such a constraint 
is external to the model defined at the quark level. Furthermore, somewhat artificially, 
the regularization in these models is applied only to the real (non-anomalous) part of 
the Euclidean quark loop term. The (finite) imaginary part is left unregularized in order 
to properly describe anomalous processes. The above-mentioned problems disappear in 
models with non-local quark interactions, which result, for example, from QCD derivations 
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of low-energy effective theories [7-9]. In such models, the quark propagators have 4- 
momentum-dependent masses which regularize the theory in the Euclidean domain. This 
regularization is in fact far more satisfying than in local models. Indeed, it is no longer 
required to regularize differently the real and imaginary parts. Furthermore, anomalous 
processes turn out not to depend on the cut-off, and are properly described [10-12]. 
Finally, the momentum- dependent regulator makes the theory finite to all orders in the 
1/Nc expansion. This is in contrast to local models, where inclusion of higher-order- loop 
effects requires extra regulators [13], and the predictive power is thereby reduced. The 
price to pay is a greater difficulty in calculating hadron properties. In view of the afore- 
mentioned advantages, it is a challenge to discover whether solitons exist in these theories. 
The purpose of this Letter is to report that indeed stable solitons do exist in a chiral quark 
model with non-local interactions. 

Our model is defined by the Euclidean action: 

/ = -Tr log {-id^-i^ + m + r$r) + ^ J d^x [S^ + P^) (1) 

where $ = 5 + i'^^Pa^a is the SU{2) chiral field. The chiral field is the dynamical variable 
of the system and it is local in coordinate space. The trace is in color, isospin, Dirac 
indices and space-time. The average u and d current quark mass is denoted by m. The 
coupling constant G has the dimension of inverse energy. The operator r is the regulator, 
chosen to be a scalar function diagonal in the 4-momentum space: {k\r\k') = 6{k — k')rk- 
It introduces a non-locality in the interaction between quarks and the chiral field. The 
particular way of introducing the non-locality is the one which results from instanton- 

induced interactions [7]. It has also been used in Ref. [14] .0 We choose the Gaussian 

_ fc^ 

form Tfc = e ^ , where A is a cut-off parameter. The form of r is not critical for results 
presented here. 

Our treatment of model (1) is the same as in other soliton calculations in effective chiral 
models: we treat the chiral 5* and Pa fields classically, which amounts to keeping the 
leading-order contribution in the number of colors, Nf.- For sufficiently large values of G the 
model leads to dynamical chiral symmetry breaking. The vacuum value of the scalar field, 
5*0, is determined by the stationary-point condition 5//5S'|^^^^ p^^g = 0, which explicitly 



^ With help of the semi-bosonization technique, it is straightforward to show that the model 
(1) is equivalent to the model of Birse et al. [15,16], defined by introducing a non-local separable 
4-quark interaction of the form 6{pi + P2 + Ps + Pi)fp^rp^rp^rp^, where pi are the (incoming) 
momenta of quarks. The meson fields of Eq. (1) are related to the quark fields of Refs. [15,16] 
in the following way: 

{S{x),Pa{x)} = -G' j e-(^'i+^^)-V^(pi)r,,{l,z75ra}r,,V'fe). 
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nves: 



The pion mass, = 139 MeV, is identified with the location of the pole of the pion 
propagator, and the pion decay constant, = 93 MeV, is easily obtained via the Gell- 
Mann-Oakes-Renner relation [15,16]. We keep the leading terms in the chiral expansion: 

Fl = 2KNjf + O p) , = ^ + O p) , (3) 



where 



dAk rl , „ r d.k r^^^, +fc 



and / = (4) 



i27rr k^ + {rlSoY J {27,^ {k^ + (rlSo) 



The model (1) has three parameters: G, A, and m. The stationary-point condition (2) and 
the equations (3) allow us to express them in terms of one parameter, Sq. Our numerical 
results are presented for various choices of 5*0.0 

The static meson fields are determined self-consistently by solving the Euler-Lagrange 
equations deduced from the action (1). In order to evaluate the trace in functional space 
it is convenient to introduce the energy- dependent Dirac Hamiltonian, h{y'^), given by 

h =-ia-V + pr - V^) $ r (z/^ - V^) + /3m, (5) 



with eigenstates and eigenvalues e\ {v"^) satisfying h (z/^) = ex (z/^) \Xu)- We obtain 

^c(;?^ p _ (A.|r|f) {(3, i^^Tg} {x\r\K) 



The integration contour over u in Eq. (6) has to be chosen in such a way as to include 
the valence quark orbit. This point requires a careful explanation. It is well known that in 
hedgehog solitons a 0+ (grand-spin zero, positive-parity) orbit becomes bound and well 



^ It should be noted that unlike the NJL model with local interactions, Sq does not have the 
interpretation of the constituent quark mass. It does not correspond to the location of the pole 
in the quark propagator, which is obtained by the condition k'^ + (r^S'o + rn)'^ = 0. In fact, in 
the vacuum sector, for sufficiently large 5*0 the quark propagator has no poles on the physical 
axis, i.e. for real negative fc^. See the discussion following of Fig. 2. 
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isolated from other orbits [17]. This feature persists in the present model. We search for 
a pole of the quark propagator in the complex v plane by solving the equation 

%v + eo (/) = 0, (7) 



where eo (z^^) is the eigenvalue of the valence 0"*' orbit calculated with the Dirac Hamilto- 
nian h (z/^). The value of u which solves Eq. (7) is denoted by iCvai, and the corresponding 
orbit by |val). We find that it is always possible to find a solution of Eq. (7) in the back- 
ground of soliton fields, even for choices of model parameters for which the quarks do 
not materialize on-shell in the vacuum [3,14,18]. The solution of Eq. (7) corresponds to 
a pole of the quark propagator which lies on the imaginary [i.e. physical) axis close to 
u = 0. It is well separated from other (unphysical) poles in the complex plane. Such poles 
are associated with the presence of the regulator [19]. In all cases described in this Letter 
Cvai > 0. The inclusion of the valence orbit in Eq. (6) yields 



{S (x) , Pa (x)} = G^(val|r|f) {f3, i(3^,Ta} (x|r|val) 



diy_ {Xu\r\x) {/?,i/?75ra} (x|r|A,,) 
27r V + ex (jy^ 

The second term is the Dirac sea contribution, where the integral over u is performed 
numerically along the real u-axis. Note that in the absence of the regulator (r = 1) 
the eigenvalues are energy- independent, the integral over u can be done via the Cauchy 
theorem, and the Dirac-sea part of Eq. (6) reduces to the usual formula expressed by the 
sum over the negative-energy spectrum. 

There is a non-trivial problem associated with theories with non-local interactions: Noether 
currents, in particular the baryon current, acquire extra contributions induced by the 
presence of the regulator r [10,16,15]. By performing the gauge transformation, ip (x) 
Qi''n{x)ijj (^x), the quark part of the action is transformed into: 

- Tr log e-*" (-29^7/. + r^r) e*" = -Tr log (^-td^^^ + + e-*V$re''' j (9) 



The extra contributions are due to the term e *'?r$re*''. When these are taken into account, 
the expression for the baryon number becomes: 



B 



1 «(„) _ 1 f,,,,^ . + ^ 



I d^T. ^ \ 'Uy (10) 



The non-local regulator contributes the extra term — ^ui which causes the residues of 
poles of the quark propagator to be normalized to 1. As a result, the baryon number is 
properly quantized, independently of the shape of the regulator.[^ 

^ Note that the baryon number in not quantized (although it is enforced to be equal to one) 
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The energy of the time independent sohton can be deduced from the euchdean action (1): 

1 T i + — 1 /• 

— oo 

where vac denotes the vacuum subtraction. The same expression can be derived from the 
Noether construction of the energy-momentum tensor. 
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Fig. 1. Self consistent fields for 


= 400 MeV 



Our numerical calculations were performed with the hedgehog-shaped chiral field, in which 

S{x) = 5'(|^|) and Pa{x) = XaP{\x\). The Euler-Lagrangc equations (8) are solved by 
iteration. The Kahana-Ripka spherical plane- wave basis [17] is well suited to calculate 
the matrix elements of the Dirac Hamiltonian h{v'^). The quark orbits arc calculated by 
diagonalizing h (z/^) in this basis, successively for each value of the energy u. We find 
convergence without forcing the fields to remain on the chiral circle, as was the case of 
solitons in Refs. [1-3]. Figure 1 shows the self-consistent fields obtained for 5*0 = 400 MeV. 
The departure from the chiral circle is substantial. Increasing 5*0 brings the fields closer to 
the chiral circle. It is important to stress that our non-local rcgularization cuts off effects 
due to high gradients in the chiral fields. Such effects were responsible for the instability 
in the case of the proper-time rcgularization [5,6]. 

Figure 2 shows that our soliton is energetically stable above Sq ~ 280MeV. We plot the 
soliton mass (dot-dashed line) and 3 times the vacuum quark mass, 3Mg (solid line). 

in the approach of Ref. [20], where the proper-time regularization is applied to the imaginary 
(anomalous) part of the effective action. 
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Fig. 2. The energy of the sohton (dot-dashed hne) and the three free-space quark mass (soUd 
Unc) plotted as functions of the parameter Sq. Above Sq ~ 280MeV the sohton is energetically 
stable. 
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Table 1 

Properties of the self-consistent soliton solutions 

The quantity Mq is defined as the pole of the quark propagator in the vacuum, i.e. 
corresponds to the solution of the equation A;^ -|- (r^-S'o + fnf\k'^=-M'^ — 0- ^^^^^ solutions 
to this equation exist only below the critical point denoted by a circle at 5*0 = 309MeV. 
Thus in the region 280MeV < 5*0 < 309MeV the mass of the soliton is lower than the 
mass of three quarks, hence the soliton is bound. Above 5*0 = 309MeV the soliton is also 
stable, since there are no asymptotic quark states (no real quark poles) into which it could 
decay to. Below -S'o = 280MeV the sohton is not energetically stable. 

Table 1 lists several soliton properties with various values of the vacuum scalar field. Pa- 
rameters A and m are determined from (3), and G from the stationarity point condition 
(2). Different contributions to the sohton energy are given: the energy of the valence orbit 

obtained from (7), the Dirac sea contribution (the second term in (11)). The contribu- 
tion of the third term in (11) is not explicitly given since it can be deduced from the 
total energy and other components given in Table 1. The soliton mean square baryon 
radius (r^) is calculated from the baryon density. Its main contribution is due to the 
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valence quarks, and the Dirac sea contributes only about 2 %. The axial-vector charge 
of the nucleon, g^, is deduced from the asymptotic behavior of the pion field, given by 
P(r) ~ —D{l/r^ + m^/r)exp(— m^r), and equals to [21,22] = SttDF^/3. For 5*0 above 
300 MeV the value of qa is somewhat smaller than the experimental value of 1.26, but 
50 % larger than the leading- A^c value obtained using proper time regularization and con- 
straining the chiral field to remain on the chiral circle [1]. Our value of qa is intermediate 
between too low Skyrme values (on the chiral circle) and too high values obtained in 
T. D. Lee solitons (where the pion field vanishes) which are more akin to the MIT bag 
model. 

Our soliton is a mean-field hedgehog solution which, of course, should not be identi- 
fied with the nucleon. Quantizing the collective coordinates, or projecting the soliton 
wavefunction on the subspace with nucleon quantum numbers, will reduce its energy by 
eliminating the spurious rotational and translational motion. It is therefore not undesir- 
able that the energies in Table 1 be higher than the experimental nucleon mass. We may 
similarly expect that the large radius of the soliton will be reduced when the center-of- 
mass corrections are calculated. The aim in this Letter was to investigate the mechanism 
of soliton formation and its stability in chiral models with non-local regulators, rather 
than to make detailed predictions for various nucleon observables. Further properties of 
the soliton and a more detailed account of the calculation will be published shortly. 
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